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Resonant Transfer of Excitons and Quantum Computation 
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Resonant energy transfer mechanisms have been observed in the sensitized luminescence of solids, 
in quantum dots and in molecular nanostructures, and they also play a central role in light harvesting 
processes in photosynthetic organisms. We demonstrate that such mechanisms, together with the 
exciton-exciton binding energy shift typical of these nanostructures, can be used to perform universal 
quantum logic. In particular, we show how to generate controlled exciton entanglement and identify 
two different regimes of quantum behaviour. 
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The Forster energy transfer was first studied in the 
context of the sensitized luminescence of solids 0, , in 
which an excited sensitizer atom can transfer its excita- 
tion to a neighbouring acceptor atom, via an interme- 
diate virtual photon. This mechanism is also responsi- 
ble for photosynthetic energy processes in antenna com- 
plexes, biosystems (BSs) that harvest sunlight 0\. More 
recently, interest has focussed on such a transfer in quan- 
tum dot (QD) nanostructures Q and within molecular 
systems (MSs) 0. In this Letter we give a general scheme 
for quantum computation that can be implemented in 
different nanostructures (NSs) by exploiting the Forster 
and exciton-exciton interactions. Thus, methods for con- 
trolled generation of exciton entanglement that use both 
diagonal and off-diagonal interactions are given. 

Consider two coupled generic NSs (scalability will be 
addressed later). We assume that the excitations of each 
NS are charge neutral (i.e., of an excitonic nature) and 
that they can be produced by optical means [§]. We 
also assume that tunnelling processes between them may 
be neglected, but that there is a strong exciton-exciton 
electromagnetic coupling. Our two-level — qubit system 
is represented in each NS by a single low-lying exciton 
(qubit state |1)) and the ground state (qubit state |0)). 
Then the interaction Hamiltonian can be written in the 
computational basis ({|00), |01), |10), |11)}, with the first 
digit referring to NS I and the second to NS II) as follows: 
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where the diagonal interaction Vxx is the direct Coulomb 
binding energy between the two excitons, one located on 
each NS, and the off-diagonal Vp denotes the Coulomb 
exchange (Forster) interaction which induces the trans- 
fer of an exciton from one NS to the other. These are 
the only Coulomb interaction terms which act between 
the qubits. ujq denotes the ground state energy and we 
define Ao = u>i — 0J2 as the difference between the exci- 
ton creation energy for NS I (u>i) and that for NS II (0^2)- 



Energy contributions due to spin singlet-triplet splittings 
do not significantly affect the present gating scheme and 
such effects are dealt with elsewhere |7J. 

The eigenenergies and eigenstates of the interacting 
qubit system are Eqq = loq, Eqi = loq + uj\ — -^(1 — A), 
E w = to +u 1 — -^-(l+A),E n = loq +uj 1 +lo 2 + Vxx; and 
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where A = + 4(Vf/A ) 2 , 
TJ/2A (w Vp/A for Vf/A < 1) and 
C2 = -J {A — 1)/2A The eigenenergies in the absence 
and presence of interqubit interactions are displayed in 
Figs.^a) and (b) respectively. Fig.^c) shows c\ and c 2 
as a function of the ratio Vf/Aq. 

Single qubit operations can be achieved by induc- 
ing Rabi oscillations in the excitonic system (e.g., see 
Ref. 11). The Vxx and Vp interactions lead to two 
regimes for achieving quantum entanglement. First, 
if the ratio Vf/Aq ^> 1, then after selectively excit- 
ing NS I and creating |10), the system will naturally 
evolve into one of the maximally entangled Bell states 
^|(|10}±i |01) ) and this evolution could then be stopped 
by suppressing the Forster coupling (see later). Second, 
if Vf/Ao <C 1, the computational basis states are essen- 
tially the eigenstates of the system, and we can exploit 
the Vxx term, together with appropriately tuned laser 
pulses, to perform two-qubit logic gates and hence to gen- 
erate controlled entanglement. From Fig.^b) we can see 
how to implement the CNOT gate; e.g., CNOTi 2 (|l) |0) h-» 
|1) |1)) can be achieved by illuminating the qubit system 
in the state 1 10) with a 7r-pulse of energy e±2- If we start 
in the ground state and first apply a tt/2 or 3ir/2 pulse at 
energy wi , we create the states -^(|00) ± 1 10) ) ; if we now 
apply a 7r £l2 pulse, we generate the maximally entangled 
states ^75 (| 00) ± |11)). 

In order to illustrate our ideas, we shall first calculate 
the strength of the interactions in a coupled QD molecule 
(in a previous study of QD excitons for quantum comput- 
ing (QC) 0, the off-diagonal interaction terms were ne- 
glected). We start with the direct Coulomb interaction, 
Vxx- Different dot geometries (e.g. spherical, pyramidal, 
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cuboidal) can be used to implement the scheme. In this 
Letter we assume that the dots are square-based cuboids 
and that the potential energy V of both electrons and 
holes increases abruptly at the cube boundaries where 
the semiconductor bandgap changes, and V = inside 
the cubes (see Fig. EJa)). This type of potential has the 
advantage of both a well defined dot size in all three di- 
mensions, and of bound and unbound solutions in each 
direction (this is in contrast with the parabolic poten- 
tial considered in Ref. [jj). The wavefunctions for single 
particles may be expressed in the envelope function ap- 
proximation as ip p (r) = (j>p(r)Up(r) |l0|, where <fi p (r) is 
an envelope function describing the changing wavefunc- 
tion amplitude of confined states for particle type p over 
the dot region, and U p (r) is the Bloch function which 
has the periodicity of the atomic lattice. In the effective 
mass approximation, the envelope functions are solutions 
of the single particle Schrodinger equation for the poten- 
tial V p , and for the effective mass m* of particle p (see 
Fig. Et a ))- These solutions may be obtained by expand- 
ing the Hamiltonian in a set of envelope basis functions 
of the form f2(r) = £ x (%)£y{y)£z(z), (the are the so- 
lutions of a one-dimensional square well potential) and 
then diagonalizing (lf)j |. 

We now combine the ground states of electrons and 
holes into a Slater determinant type wavefunction which 
represents one exciton on each dot: 

*xx = ^(n)^h(ra)V?('a - R)V£(r4 - R)] , (2) 

where A indicates that the wavefunction has overall anti- 
symmetry, R is the vector connecting the two dot centres, 
ri and r 3 represent the position vectors of electrons rela- 
tive to the centres of dot I and dot II respectively and r 2 
and r4 are the equivalent vectors for holes. The expan- 
sion of the associated Coulomb operator Vxx in a Taylor 
series about R to lowest non-zero order yields 

Vxx = ^5 | Pi ■ Pn - ^(Pi ■ R)(Pn • R)| j (3) 
where k = -r—, £ r denotes the dielectric constant of 
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the medium (e r = 10 throughout our discussion), pi = 
e(ri — 1*2), and pn = e(r3 — 1*4) are the overall dipole 
moments on dot I and II respectively. To evaluate the 
matrix element (^xxl^xxl^xx), Pi and pn in Eq. © 
are replaced by their expectation values for the wave- 
function, Eq. J5J. This procedure gives rise to a direct 
term and exchange terms. The exchange terms are zero 
in the absence of wavefunction overlap between dots. The 
direct term is obtained through the use of the envelope 
function approximation which leads to the expectation 
value (ri) = J </>e*(ri)ri0* (ri)dri, and similar expres- 
sions hold for the other position expectation values. For 
a cuboidal dot, where the electron and hole wavefunc- 
tions have a definite parity about the dot centre, Eq. (jjJJ 
implies that the exciton-exciton coupling is zero. How- 
ever, when this symmetry is broken, the electron and hole 



localize in different parts of the dot and the dipole mo- 
ment is non zero. This occurs, for instance, i n pyra midal 
dots 0| or when an electric field is present |flL Il2j. 

We have simulated the effect of applying an electric 
field (along x) in our model and the results are displayed 
in Fig.[2Ib), where the exciton dipole moment pi is plot- 
ted as a function of the dot size and the field strength 
E. Two cases are considered: the first is that of a cubic 
dot (a = h/2) and the second is that of a flat cuboidal 
dot (a = 5/i). Vxx is obtained by using the size of pi for 
each dot, and Fig. EJc) shows Vxx (normalized by R~ 3 ) 
calculated for both geometries, with a = b. It is interest- 
ing that the size of the induced dipole depends strongly 
only on the length, a, of the dot in the direction of the 
applied field, and at large field is given approximately by 
Pi ea; this limit should be valid for any NS (such as the 
BSs and MSs described later) with a well defined length- 
scale, and not just for QDs. The cuboidal structure is 
more typical of Stranski-Krastanow self assembled dots, 
where typically R = 5 nm, a = 10 nm, and b = 8 nm, 
and hi = hi = 2 nm. In a field of 100 kV/cm, these 
parameters give Vxx = 120 meV, which would result in 
a lower time limit for the gate operation of around 10 fs. 
This is relatively short; decoherence times on the order of 
nanoseconds have been observed for uncoupled dots . 

Let us now consider the Forster (off-diagonal) cou- 
pling Vp in QDs, which may be expressed as a matrix 
element of the Coulomb operator between excitons lo- 
cated on each of the two dots. By Taylor-expanding this 
expression around R to lowest non-zero order |2j, we 
obtain V F = ^3 «n) ■ (rn) - Jr«n> • R)((r„) • R)). 
Here e(rj) (i = I, II) represents the matrix element of the 
position operator between an electron and a hole state on 
dot i. The expression is therefore equivalent to the in- 
teraction of two point dipoles, one situated on each dot. 
However, its magnitude is quite different from Vxx- By 
again employing the envelope function approximation, we 
can rewrite the above equation as 

Vf = ^JojOn ((if) • (rft) - A«rf) • R)((rf : ) • R)) , 

where (rf) = j „ J7g(r)r[/^(r)dr is the interband posi- 
tion matrix element for the atomic part of the wavefunc- 
tion for dot i, and O, = J (j) l c (r)(j) l h (r)dr is the overlap of 
electron and hole envelope functions on dot i. Thus, the 
effects of the QD size and shape (which determine the 
overlap integrals) are separated from the effects of the 
material composition of the dot (which determine the in- 
terband position operator). For other, smaller, NSs, the 
envelope function approximation does not apply. How- 
ever, in these cases the above analysis is still valid if the 
overlap integrals are set to unity, and if the NS's entire 
wavefunction is represented by the atomic basis part of 
the general wavefunctions described above. 

We now introduce a simple Kronig-Penney model to 
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describe the atomic basis, and we assume that the infi- 
nite Kronig-Penney quantum boxes have a width of 2x. 
This gives (rf) — 32x/9ir 2 , or Vp « i 2 , and so we plot 
Vp/x 2 as a function of R in Fig.[3Ja). The simple Kronig- 
Penney model shows how the size of the Forster transfer 
depends upon the physical size of the atomic part of the 
wavefunction. (r a ) is a widely measured quantity since it 
determines the strength of dipole allowed transitions in 
optical spectra. In CdSe QDs it can be in the range of 

An 
|4J , in atomic systems it can also be several 

and in BSs and MSs has recently been observed to 
be about 1.7 eA $\. The solid line in Fig.^a) represents 
the case when the overlap integrals are set to unity; the 
symbols represent simulations which take higher order 
terms into account and have a = b = 2 nm .7]. Fig.^b) 
shows how Vp varies as a function of dot size and confine- 
ment potential through the Ot 's of Eq. . The envelope 
functions are as described in the previous section, except 
that no electric field term has been included in the en- 
velope Hamiltonian. Thus the values of Vf should be 
regarded as an upper limit for the cuboid dot model: ap- 
plying a field will serve to reduce Vf since it decreases 
the Oi's; hence it could be used to tune Vf- The overlap 
is enhanced when there is a larger confinement potential 
and for larger dots, since in these cases the shape of the 
wavefunction is less sensitive to the effective mass of the 
particle. As an illustration of the use of these curves, let 
us assume that we have a dot system in which, as before, 
R = 5 nm, a = 10 nm, 6 = 8 nm, and h\ = hi = 2 nm. 
Furthermore, let us take the measured dipole value for 
CdSe dots of 0.9 to 5.2 eA 0. In this case, the Forster 
strength is between 0.02 and 0.6 meV, which if Ao = 
would correspond to an on resonance energy transfer rate 
of between 206.8 and 6.9 ps. This is short enough to be 
useful for QIP: decoherence times as long as a few ns [lflj 
have been observed in QDs. In MSs or BSs, the interact- 
ing units can be as close together as 1 nm; using this and 
taking a typical molecular or biomolecular dipole value 
of about 1.7 eA 0,0, we obtain an interaction strength 
of 8.3 meV (or a transfer time of 497 fs). Furthermore, 
Vp must certainly be controlled if the alternative scheme 
using Vxx is to be implemented (and therefore cannot be 
neglected as in Ref. |||). 

By also using the model outlined above to calculate A 
in QDs (shown in Fig. Ufa) &s a function of dot size ra- 
tio for quantum cubes) we can compute the size of the C\ 
component of the j^io) and | M/ oi ) states. This is shown in 
Fig-Efb), where the dependence on atomic dipole and dot 
separation are incorporated by multiplying c\ by R?x~ 2 . 
It can be seen there that a range of c\ values can be 
obtained by choosing dots with appropriate values of x, 
R and a/b. Dots with large x (> 1 nm say), small R 
(< 3 nm say) and a/b ~ 1 give a larger c\, and it is 
then more appropriate to use the Forster interaction it- 
self to create entangled states. On the other hand, dots 
with smaller x, larger R or a large mismatch in dot size 



would be more suited to the scheme which uses the Vxx 
for QC and entanglement generation. The fidelity of a 
typical operation (e.g. |11) i-> 1 10)) in this case is equal 
to 1 — c 2 — and so one must be careful when using the 
biexciton scheme to use the available parameter space 
and make sure that the Forster transfer is suppressed to 
the desired accuracy. There are other sources of decoher- 
ence in this case (e. g. t he interaction with optical and 
acoustic phonons |3,ll3j) which will reduce the value of 
the fidelity to below 1 — c\ . To minimize the effects of 
such decoherence channels, it is important to maximize 
Vxx, since this leads to an improved transition discrimi- 
nation and so to a faster gating time. This can be done 
by applying an electric field and choosing an appropriate 
dot shape, size and separation (as described earlier). It is 
then necessary to minimize the basis state mixing for the 
chosen parameters by selecting a suitable dot size ratio 
and material composition. 

Single shot qubit state measurement in QDs could be 
performed by using resonant fluorescent shelving tech- 
niques [l5j. The QD state measurement can also be 
achieved by means of projecting onto the computa- 
tional basis and measuring the final register state by ex- 
ploiting ultrafast near-field optical spectroscopy and mi- 
croscopy 0: these allow one to address, to excite and to 
probe the QD excitonic states with spectral and spatial 
selectivity. In addition, the qubit register density matrix 
can be reconstructed by measuring the QD photon cor- 
relations via standard quantum state tomography tech- 
niques [l^. Scalability of the scheme given here could 
also be possible by adopting a globally addressed qubit 
strategy |ljj on a stack of self-organized QDs |18j . 

Light-harvesting antenna complexes [3j or arrays of 
strongly interacting individual molecules 0] could pro- 
vide an excellent system in which the Forster interac- 
tion could be used for QIP tasks. They are generally 
very uniform structures, and we may compare them to 
QDs by setting a/b ~ 1 5 or Vp/A ^> 1. Then the one- 
exciton eigenstates of a two qubit system with a Forster 
coupling naturally allows the generation of the states 
^7j(|01) ± 1 10) ) , which, apart from their applications to 
quantum protocols, can be particularly useful in the fight 
against decoherence. Spectroscopic, line-narrowing tech- 
niques (e.g., hole burning and site-selective fluorescence), 
infrared and Raman experimental studies reveal that the 
main decoherence mechanisms in the antenna complexes 
arise from energetic disorder, electron-phonon coupling, 
and temperature effects Q. In this scenario, the excita- 
tions couple to an environment that typically possesses a 
much larger coherence length than the biomolecular units 
(BChl's) spacing. For example, the BChl's in the antenna 
complex LH2, which we consider as potential qubits, are 
spaced by as little as 1 nm, and hence so-called collective 
decoherence is expected to apply. In this case, provided 
that the logical qubit encoding = |01) jfc , |t), i = |10) ife 
that uses two physical (exciton) qubits can be realized 
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in the BChl's system, arbitrary superpositions of logi- 
cal qubits such as (a< IJ.^ + A \1) t )® N , i = 1,---,N, 
cti, Pi € C, are immune to dephasing noise (described 
by a a z operator |l9|). and single qubit manipulations 
can be carried out on the timescale of the Forster cou- 
pling (which as we have seen can be as short as 497 fs). 
Two-qubit logic gates can also be implemented within a 
decoherence-free subspace by using the above encoding, 
thus completing a universal set of gates 0. Initializa- 
tion of the system requires the pairing of the physical 
qubits to the logical 'ground' state \i)f N , and readout 
is to be accomplished by identifying on which of the two 
structures the exciton is. Furthermore, rings of BChl's 
appear side by side in naturally occuring antenna com- 
plexes and also display energy selectivity — smaller rings 
tend to have higher energy transitions [3j. Thus, follow- 
ing a scheme as above, it may be possible to scale up such 
biological qubits in a natural way and construct a robust 
energy selective scheme for quantum computation. 

In conclusion, we have provided a general scheme for 
quantum computation and quantum entanglement gen- 
eration that can be implemented in different NSs by ex- 
ploiting the Forster and exciton-exciton interactions. In 
particular, we have shown how such interactions can be 
manipulated in molecular, biomolecular and QD nanos- 
tructures in order to produce an accurate degree of con- 
trol for quantum logic. 
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FIG. 1: A schematic diagram of the proposed quantum logic 
gate scenario: energy levels in (a) the absence and (b) the 
presence of the interactions Vxx and Vf for nanostructures 
of different excitation frequency. £12 = u>2 + Vxx — S, 621 = 
t*ti + Vxx + 3, 5 = Vp/Ao; (c) dependence of the eigenstate 
coefficients d as a function of Vf/Aq. The inset shows the 
eigenenergies E01, and E10 of plot (b) for oji/Aq = 20. 



FIG. 2: (a) Schematic diagram of the cuboidal dot model. 
The potential inside the cuboids is set to zero, that outside is 
determined by the band offsets of the conduction and valence 
bands within the heterostructure. (b) Exciton dipole moment 
as a function of the dot size and applied electric field for 
two dot shapes. The dot parameters m e = 0.06, mj, = 0.6, 
V R = Vh = 500 meV. (c) Exciton-exciton binding energy for 
a = b and sequence of dot shapes, size and field strength as 
in (b). 



FIG. 3: (a) Dependence of the Forster interaction strength 
on the interdot separation. The solid line represents the case 
where Oi = 1 in the dipole-dipole approximation and the 
atomic dipole operator is given by 32x/9ir 2 . The symbols 
result from a full numerical simulation for a — b — 2 nm, 
and the dotted lines are the dipole-dipole predictions in these 
cases, (b) The Oi factor appearing in Eq. Q as a function 
of dot size and confinement potential. m e = 0.06 mo, ran = 
0.6 mo. A lower cut-off occurs when the ground state of the 
dot is no longer a bound state. 



FIG. 4: (a) Energy splitting Ao of the qubit exciton states 
l^oi), and 10) in the absence of the Forster interaction as 
a function of the different dots sizes a and b. The splitting is 
independent of interdot distance, (b) The size of the mixing 
component of the wavefunction, ci, as a function of the dot 
size ratio, ci has been scaled by its dependence on the interdot 
distance, R, and typical atomic spacing, x. 
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